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We consider the phenomenon of quantum charge pumping of electrons across a superconducting 
double barrier structure in graphene in the adiabatic limit. In this geometry, quantum charge 
pumping can be achieved by modulating the amplitudes (Ai and A2) of the gaps associated with 
the two superconducting strips. We show that the superconducting gaps give rise to a transmission 
resonance in the A1-A2 plane, resulting in a large value of pumped charge, when the pumping 
contour encloses the resonance. This is in sharp contrast to the case of charge pumping in a normal 
double barrier structure in graphene, where the pumped charge is very small, due to the phenomenon 
of Klein tunneling. We analyse the behaviour of the pumped charge through the superconducting 
double barrier geometry as a function of the pumping strength and the phase difference between 
the two pumping parameters, for various angles of the incident electron. 

PACS numbers: 73.23.-b,72.80.Vp,74.45.+c 



I. INTRODUCTION 

The phenomenon of quantum charge pumping corre- 
sponds to a net flow of DC current between different 
electron reservoirs (at zero bias) connected via a quantum 
system whose parameters are periodically modulated in 
tim(P2l. The zero-bias current is obtained in response 
to the time variation of the parameters of the quantum 
system, which explicitly break time-reversal symmetry. 
It is necessary to break time-reversal symmetry in order 
to get net pumped charge, but it is not a sufficient condi- 
tion. For obtaining a net pumped charge, parity or spa- 
tial symmetry must also be broken. Within a scattering 
approach, if the time period of modulation of the scat- 
tering system parameters is much larger than the time 
the particle spends inside the scattering region (dwell 
time), the adiabatic limit is reached. In this limit, the 
pumped charge in a unit cycle becomes independent of 
the pumping frequency. This is referred to as "adiabatic 
charge pumping""^. In recent years, quantum charge and 
spin pumping through various mesoscopic samples, in- 
volving quantum dots and quantum wires, have attracted 
increasing interest both theoretically "^'^^ and experimen- 
talljii2H2il^ both in the adiabatic regime and otherwise. 

The discovery of graphene, a two dimensional single 
layer of graphite, by K. S. Novoselov et al^ a few years 
ago, has led to an upsurge in the study of its traiisport 
properties, both theoretically and experimentally^MSSl^ 
The low energy quasiparticle excitations in graphene be- 
have like massless relativistic Dirac fermions. This pro- 
vides us with an experimental test bed for observing 
many well-known phenomena in relativistic quantum me- 
chanics, such as the Klein paradojP^ at low energies. In 
the recent past, a graphene-based quantum pump has 
been considered in literature I^HSU -^^i^gi-g pumped charge 
is obtained in an adiabatic quantum pump device based 
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FIG. 1: (Color online) Cartoon of the SDB structure where 
a graphene sheet is connected to two reservoirs labeled 1 and 
2. Superconducting material is deposited on top of the two 
patches labeled S and connected to contacts labeled 3 and 
4. The schematic of the potential profile seen by an incident 
electron is shown below. 



on a graphene monolayer modulated by two oscillating 
gate potentials. However, the pumped charge obtained 
in these kind of devices is quite small. 

Quantum charge pumps, using a variety of setups in- 
volving superconductors, have also been of major inter- 
est in recent years^-"^. Also, very recently, adiabatic 
charge pumping in graphene with superconductors has 
also been considerecpD. However an adiabatic quantum 
pump device based on superconducting double barrier 
(SDB) structures in graphen^^ has not yet been con- 
sidered in the literature. Motivated by this fact, in this 
article we consider quantum pumping of electrons (in the 
adiabatic limit) across an SDB structure in graphene, as 
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depicted in Fig. [T] Till date, no experiment has been 
carried out in the context of charge pumping with super- 
conducting barriers. Experimentally it might be possi- 
ble to design an SDB structure by depositing thin strips 
of superconducting material on top of a single ballistic 
2D graphene sheet at two places. This can induce a fi- 
nite superconducting gap in the barrier regions of the 
graphene sample due to the proximity of the supercon- 
ducting strips. In principle, we can explore two scenarios 
to achieve significant amount of pumped charge — (a) hy 
periodic modulation of amplitudes Ai and A2 of the gaps 
at the two superconducting barriers (SB) or alternatively, 
(b) by periodic modulation of the order parameter phases 
and (p2 associated with the two barriers. In this paper, 
we explore the first alternative, since it has been seen in 
earlier worlpP that the second alternative leads to less 
pumped charge. For free electrons in the graphene sheet, 
we show that in the Ai — A2 plane, there is a sharp res- 
onance point in the transmission probability across the 
SDB structure. This is in sharp contrast to the case of 
charge pumping in a normal DB structure in graphene 
where the transmission probability across the DB struc- 
ture does not have any resonance structure due to the 
phenomenon of Klein tunneling. In the SDB geometry, 
when we consider Ai and A2 as the pumping param- 
eters, we can always choose a pumping contour which 
completely encloses the transmission resonance. Hence 
it is possible for the pumped charge to be large if the 
resonance is sharp enough. 

This paper is organized as follows. In Sec. |llj we 
discuss the theoretical modelling of the SDB structure 
in graphene and describe the resonance structure in 
the pumping parameter Ai — A2 plane. In Sec. Ill 



we 



discuss the numerical results for the pumped charge as a 
function of the various parameters in the theory and also 
contrast it with that of normal DB geometry. Finally in 



Sec. IV we present our summary and conclusions. 



electron) can be expressed as 

V{x) = Ae^'^e(x)e(-a; + a) + Ae'"^ 

e[x - {a + L)]e[-x + {2a + L)] (1) 

where a is the width of the superconducting barrier and L 
is the distance between the two barriers. Here we assume 
that the spatial variation of the potential steps is slow 
on the scale of the lattice spacing so that inter-valley 
scattering is suppressed. 8 is the Heaviside 8-function, 
and we have taken (j>i = 4'2 = <i>, since we will not be 
looking at supercurrents (Josephson effect) in this work. 

Following Ref . |12J we use a four dimensional version of 
the Dirac-Boguliobov-de Gennes equation (DBDG^P'for 
electrons and holes which is given by 
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where, U = U{r)+Ep, and the energy e is measured from 
the Fermi level of the superconductor. We assume that 
U{r) — in the normal graphene region and C/(r) Uq, 
a constant, independent of r in the proximity induced 
superconducting region. Note that we have defined di- 
mensionless variables 
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to replace the original ones. 

The solution of the DBDG equationJ^, describing elec- 
trons and holes with incident energy e inside the normal 
graphene regions (A(i) = 0), can be written as 
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II. RESONANT TRANSMISSION IN 
GRAPHENE SDB STRUCTURE 



Quantum transport in the SDB structure in graphene 
was considered recently in Ref. W2!. Here we consider the 
same set-up, but instead of applying a bias across termi- 
nals 1 and 2, we shall use terminals 3 and 4 to change 
the amplitudes of the superconducting barriers in a time- 
dependent and out of phase fashion. The SDB structure 
is formed by depositing thin strips of superconducting 
material on top of the graphene sheet at two places. This 
induces a finite superconducting gap (Aic"^') in the bar- 
rier regions as a result of the proximity effect, (i refers 
to the index of the strips). The geometry is shown in 
Fig. [1] The spatial dependence of the order parameter 
(which also acts as a scattering potential for the incident 



h± 



cos a' 






±1072 



(5) 



where a sin "'"[fc.y/(e + 1)], a' = sin "'"[fcy/(e — 1)], fc 



;2 — fc2 and fc' 



fc2. q; is the angle of incidence 



of the incoming electron (with wave-vector {k,ky)) and 
a' is the angle of reflection of the Andreev reflected hole 
(with wave- vector {k' ,ky)). For retro Andreev reflection 
(AR), a',fc' have opposite signs from a,k whereas for 
specular Andreev reflection (SAR), they have the same 
signs. The change from retro AR (e < 1) to SAR (e > 1) 
occurs at e = 1 (in our dimensionless units). 

Similarly for the superconducting barrier regions, the 
four component spinor solutions (u, v) contain electron 
wave-functions u of one valley and hole wave-functions 
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FIG. 2: (Color online) Contours of the transmission probability \tc\^ in the Ai — A2 plane for three different values of ky. 
The steps between the maxima and minima of |tcp range approximately 1.0 to 0.1, 0.6 to 0.05 and 0.3 to 0.01 for the 
three cases ky = 0.0,0.3 and 0.75 respectively. In (a), (b) and (c) a/L = 0.017, Ao = 40.0, f/o = 10.0,01 = <t>2 = 0.0 and 
e = 0.00142,0.00121,0.00082 for the three contour plots respectively. The black circle represents the pumping contour for the 
parameter values uj — 1.0, P = 30.0 and 77 — 47r/9. 



V of the other valley. The DBDG equation can now 
be solved for any arbitrary energy e and the four solu- 
tions inside the superconducting barriers are given in the 
preprint version of Ref . i43j 
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where the subscripts 1/2 refers to the upper and lower 
signs on the RHS respectively, and similarly for 3/4 and 
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Here, we have not taken the limit U ^ A, e. We have 
also obtained the solution for both right- moving and left- 
moving electrons and holes. 

Solving the BdG equation in the normal and super- 
conducting regions, we obtain the net quantum mechan- 
ical amplitudes for reflection, transmission (co-tunneling 
(CT)), AR (and SAR) and crossed Andreev reflection 



(CAR) (and specular crossed Andreev reflection (SCAR)) 
of an electron incident on the SDB structure, after it has 
traversed both the barriers. 

In the non-relativistic case, a DB structure always lead 
to resonances and this affects the transmissions and the 
reflections through the system. For relativistic electrons, 
the standard paradigm is that one c annot obtain con- 
fined carrier states for normal incidenc^^HS! (j^g Klein 
tunneling. However, for relativistic electrons with super- 
conducting barriers, even for normal incidence, discrete 
Andreev bound levels can clearly lead to resonant trans- 
missions in a SDB structure in graphen^^H 

Here we use the standard wave-function matching tech- 
nique to solve such scattering problems to obtain all the 
four quantum mechanical amplitudes across the SDB ge- 
ometry in graphene. To the left of the SDB structure, 
with an incident electron from the left, the wave-function 
can be written as 



(10) 



and to the right of the SDB structure, the wave-function 
can be written as 



h+ 



(11) 



Hence matching the wavefunctions in the normal and 
proximity induced superconducting regions (Eq.(H][7|) at 
the four normal— superconductor (NS) interfaces (ir — 
0,a,a-\- L, 2a -\- L) forming the SDB structure, we obtain 
sixteen linear equations. Numerically solving these six- 
teen equations, we obtain the four amplitudes Tc, tac, tc 
and tAc, for the SDB structure, for an incident electron 
with energy e below the gap A^"^ . Here Tc is the reflection 
amplitude, tac is the AR (and SAR) amplitude, is the 
CT amplitude and tAc is the CAR (and SCAR) amplitude. 
Note that we distinguish between electron and hole pa- 
rameters, and hence the four amplitudes will be different 
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for incident electrons and holes. Moreover, in our numer- 
ical analysis we do not distinguish between specular and 
retro Andreev reflections. 



III. NUMERICAL RESULTS FOR THE 
PUMPED CHARGE 




FIG. 3: (Color online) The value of the pumped charge 
Q in units of electron charge e, for pumping in Ai — A2 
plane, is shown as a function of the pumping strength P 
for three different values of ky. Here a/L = 0.0f7, cj = 
1.0,77 = 47r/9,Ao = 40.0, Uo = 10.0, <^i = <^2 = 0.0 and 
e = 0.00142, 0.00121, 0.00082 for the three plots respectively. 

In our numerical analysis for calculating the pumped 
charge through the SDB structure, we choose the ampli- 
tudes of the order parameters of the two superconducting 
strips (Ai and A2) as our pumping parameters. Here, Ai 
and A2 are taken to oscillate with the frequency to, with 
a modulation parameter P, and a phase difference 2r] 
between them, i.e., we choose 



Ai = [Ao + P cos{Lut + T])]e 
A2 = [Ao + P cos{ujt - Tj)]e 



(12) 



Aqe is the mean value of the amplitude around which the 
two pumping parameters are modulated with time and 
P is called the pumping strength. Here e is the energy 
of the incident electron. It is adjusted to be close to an 
Andreev bound level, so that the SDB is close to a res- 
onance. Hence, effectively by varying P, for fixed e we 
vary the ratio e/A. By fixing e to be close to the reso- 
nance, we maximise the pumped charge. The presence of 
two time-varying potentials with a phase difference be- 
tween them explicitly violates parity, which is a necessary 
condition for obtaining pumped charge. The frequency 
of the potential modulation is kept small in comparison 
to the characteristic times for traversal and reflection, so 
that the pump is in the adiabatic limit. 

In Fig.l2[a), we plot the transmission probability (|<cP) 
in the Ai — A2 plane for various values of the transverse 
electron momentum ky. Unlike the case of normal DB , 
where there is no resonance because of the phenomenon 
of Klein tunneling, for the SDB , |icP = 1 resonance does 



occur for normal incidence of electrons, i.e. ky = 0. These 
bound states are produced due to the superposition of 
both electron and hole states and not just from any one 
of them. Hence these resonances are also different from 
the resonances that occur in normal double barriers in 
other materials. For nonzero values of ky, |icP — 1 gets 
damped as normal reflection and CAR (and SCAR) also 
take part in transport with nonzero values^ as shown in 
Figs.[2];b) and (c). 




FIG. 4: (Color online) The value of the pumped charge Q 
in units of electron charge e, for pumping in Ai — A2 plane, 
is shown as a function of the phase difference between the 
two pumping parameters 77 for three different values of ky. 
Here a/L = 0.017, a; = 1.0, P = 30.0, Aq = 40.0, Uo = 10.0, 
= 02 = 0.0 and e = 0.00142,0.00121,0.00082 for the three 
plots respectively. 



U sing th e modified version of the Brouwer's for- 
mul£pI2SEni^ the pumped charge through the graphene 
SDB structure can directly be obtained from the para- 
metric derivatives of the 5-matrix elements. For a single 
channel iS-matrix, we have 



{\rAc\^^ + \tAcm 



cos a 



(13) 



where 9, x, /3 and 7 correspond to the phases of the 
reflection, transmission, AR (and SAR) and CAR (and 
SCAR) amplitudes respectively. Note the negative sign 
in Eq |13[ which results from the fact that vac aud tAc 
correspond to the conversion of an electron into a hole. 
Thus, the pumped charge through the SDB structure in 
graphene is directly related to the amplitudes and phases 
that appear in the iS-matrix. Inserting the unitarity re- 
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1 in Eq 13 we finally 
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FIG. 5: (Color online) The maximum value of the pumped 
charge Qmax in units of electron charge e through the 
SDB structure, for pumping in Ai — A2 plane, is shown as a 
function of the incident angle a. 



obtain 



Q = 



2^ Jo 



dt 



\rc^ \ 9 cosa + ficosa 
+ l^cP (^X cos a + /? cos a' 
+ I^AcP — 7^ cos a' ~ P cos a' 



(14) 



Eq |14| is the working formula for pumped charge in our 
case. Note that if we substitute ky — 0, (i.e. a = a' = Q) 
in Eq{T4j then the formula reduces to 



e 

2^ 



dt 



\tc\\x + P)-P 



(15) 



which is precisely the modified Brouwer's formula used 
in Ref .40, for the quantum wire case. In Eq |15[ the first 
term is called the dissipative part" and the second term 
is known as the "topological part" and depends entirely 
on the time derivative of the AR phase. 

The pumped charge is obtained by using Eq(T4]with Ai 
and A2 as the two pumping parameters. These two pa- 
rameters are varied by periodically varying the top gate 
voltage which controls the Fermi energy of the electrons 
in the superconducting region. Thus essentially, as men- 
tioned earlier, e/A is varied for the two barriers peri- 
odically. Using EqfTil we obtain the pumped charge for 
various parameters of the system. In Fig.[3]and Fig.|4j we 
show the behaviour of the pumped charge as a function of 
the pumping strength P and the phase difference between 
the two pumping parameters 77 respectively, for three val- 
ues of ky. Note that in Fig. |3][a), the pumped charge 
increases with the increase in the pumping strength P 
as a larger value of P corresponds to a larger pumping 
contour which encloses more and more of the resonance. 



Note also that the |tcP = 1 resonance at fcj, = is not 
a sharp resonance; it has a finite width because it also 
has contribution from the " dissipative parf as shown in 
Eq |15[ and is not purely topological. The " dissipative 
parf effectively reduces the pumped charge from integer 
values but we still obtain a fairly large value of pumped 
charge - where by large we mean that the pumped charge 
is a sufficiently large fraction of unity (roughly between 
0.1 and 0.35). This is in sharp contrast to normal DB in 
graphene where the pumped charg^^ZI jg very small (in the 
range of 10""*, because there is no resonance due to Klein 
tunneling. However for oblique incidence i.e. ky ^ Q (see 
Figs. [3][b) and (c)), we obtain relatively smaller values 
of pumped charge as in this case normal reflection, CT, 
AR (and SAR) and CAR (SCAR) also contribute to Eq[l4 
and the interplay between all the quantum mechanical 
amplitudes and their phases result in smaller value of 
pumped charge. 

In Fig. [4] we show the oscillatory behaviour of pumped 
charge as a function of the phase difference ry between 
the two time varying parameters. Here also we note 
that we obtain smaller values of pumped charge through 
the SDB geometry as we vary ky from normal incidence 
to oblique incidence. Note also that in Fig. |4j for all 
three values of ky, the pumped charge becomes maxi- 
mum around 77 ~ ±7r/^. 

Finally, we also show the systematic behaviour of the 
maximum value of pumped charge through the SDB ge- 
ometry as a function of the incident angle a of the in- 
cident electron in Fig. [5] This behaviour clearly shows 
that the maximum value of the pumped charge becomes 
smaller as we vary a from normal incidence [a — 0) to 
oblique incidence (a 7^ 0). 



IV. 



SUMMARY AND DISCUSSIONS 



To summarize, in this paper, we have studied adiabatic 
quantum pumping through an SDB structure in graphene 
and have shown that in the Ai — A2 plane, pumped 
charge can be large (around 0.2- 0.3) in magnitude for 
normal incidence. This is in contrast to normal double 
barriers in graphene, where the pumped charge is small 
(around 10^**) due to the phenomenon of Klein tunneling. 
We also show that in the Ai — A2 plane, transmission res- 
onances (|icP = 1) occur due to the formation of Andreev 
bound states between the two superconducting barriers. 
When the pumping contour encloses this j^cP = 1 res- 
onance, we obtain significant amount of pumped charge 
for normal incidence {ky = 0). However, the transmission 
resonances get damped as the incident angle increases 
due to reflection and SAR. Hence, we obtain much smaller 
values of pumped charge for oblique incidence. We have 
also studied the pumped charge through the SDB struc- 
ture as a function of the pumping strength (P) and the 
phase difference (77) between the two pumping param- 
eters. The most interesting study is the evolution of 
the maximum value of pumped charge through our pro- 
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posed geometry as a function of the angle of incidence of 
the incoming electron. This shows monotonic decrease of 
pumped charge as we increase the angle of incidence of 
the incoming electron. 

Similar behaviour has also been predicted for many 
other systems where one studies quantum pumping 
through nanostructures. Integer pumped charge has been 
shown for pumping thro ugh quant um dots^^^^^ as well as 
through Luttinger liquidsm^ESIini in more recent times, 
similar behaviour of pumped charge has been predicted 
in graphene NIS junctions'*^ and in an InAs Josephson 
pumplSl. 

As far as the practical realization of such an SDB struc- 
ture in graphene is concerned, it should be possible to 
fabricate such a geometry by depositing thin strips of 
a spin singlet superconductors (like Al or Nb) on top 
of a graphene sheet!^ at two places. The width of the 



strips should be of the order of the superconducting 
phase coherence length (10 — 15nm in case of Nb) for 
CT and CAR (and SCAR) to take place. In our geom- 
etry, pumped charge can be obtained by periodically 
varying the top gate voltage which controls the Fermi 
energy of the electrons in the superconducting region 
which amounts to varying e/A for the two barriers peri- 
odically. The pumped current should be in the range of 
pico-amperes when the pumping frequency is the order 
of a few MHz, and should be experimentally measurable. 
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